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We consider the double trapping reaction A + _B—>_B,_B + C^C in one dimension. The survival 
probability of a given A particle is calculated under various conditions on the diffusion constants of 
the reactants, and on the ratio of initial B and C particle densities. The results are shown to be of 
more general form than those obtained in previous work on the problem. 
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I. INTRODUCTION 

The simple trapping reaction A + B ^ B has been 
studied intensely over the last two decades, motivated 
in part by studies of the related two-species annihilation 
problem, A -\- B 0, introduced in the classic paper 
of Toussaint and Wilczek In one spatial dimension 
the survival probability, QA{t)i of a single diffusing A- 
particle moving in an intially homogeneous background 
of diffusing B particles has been proved rigorously to de- 
cay as QA{t) ~ exp(— Ai-^/^) 0. Only recently, however, 
has an exact expression for the coefficient A been derived 
0: A = pbVDb, where pB and Db are the den- 

sity and diffusion constant respectively of the B parti- 
cles. The result has been generalized pj to all (continu- 
ous) dimensions d < 2. A curious feature of the result is 
the absence of any dependence on the diffusion constant. 
Da, of the A particle (at least in the leading asymptotics 
- Da does appear, however, in subdominant terms j^). 
For d > 2 the rate-equation approach, which predicts a 
simple exponential decay, is qualitatively correct 0]- 

In contrast to the simple trapping reaction, more com- 
plex trapping sequences have received little attention. 
This paper is devoted to the double trapping reaction 
A + B^B,B + C'^C. We try to compute the frac- 
tion of A particles remaining at time t, or equivalcntly, 
the probability, QA(t), that a single A particle has sur- 
vived up to time t, given that it is initially surrounded by 
a sea of Poisson distributed traps, the B particles. The B 
particles themselves are initially surrounded by Poisson 
distributed traps - the C particles - and thus disappear 
from the system at each time step with a probability to 
be determined. By "Poisson distributed" we mean that 
the probability to find N traps in an interval of length L 
is [{pB,cL)'^ /N\] exp(— pB.ci), where ps.c is the density 
of S or C particles. Since the C particles diffuse inde- 
pendently, and do not disappear from the system, they 
remain Poisson distributed at all times. The spatial dis- 
tribution of the B particles may change, however, as B 
particles are absorbed. 

Of particular interest is the limit t oo. Since the B- 
particle density ultimately decays to zero, the function 
QA(t) will approach a non-zero limit Qa^oo) which, on 
dimensional grounds, can only depend on the density ra- 
tio pbq/pc, where pbq = Ps(0) is the initial i3-particle 



density, and on the ratios Da/Db, Db/Dq of diffusion 
constants. The calculation of Qa(oo) is our main goal. 

This double trapping problem has been studied using 
a mean-field (i.e. rate equation) approach, which should 
again be qualitatively valid for d > 2, and by a version 
of the Galanin model for d = I 0] . The latter approach 
predicts the following limiting result: 
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Pbq 
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Db+Dc 



(1) 



This paper considers various limiting regimes of the gen- 
eral problem in which exact results or bounds can ob- 
tained for comparison with Eq. iQJ . We also obtain some 
results on the asymptotic form of the time-dependence of 
QA{t), i.e. on the manner in which the asymptotic limit 
is approached. 



II. ANALYSIS OF THE SURVIVAL 
PROBABILITIES 



Each particle i diffuses according to the Langevin equa- 
tion Xi — r]i{t), where rii(t) is Gaussian white noise with 
mean zero and correlator {rii{t)Vj{t')) = 2D6ijS{t — t'), 
where D = Da, Db, Dq are the diffusion constants 
for the three reactant species. From this one can de- 
rive in standard fashion a backward Fokker-Planck equa- 
tion governing the time evolution of the particle's survival 
probability Q{xi,t): 



^Q{x^,t) ^^fQ{x^,t) 



dt 



(2) 



where Xi is the position of the ith particle at time t — Q. 

We consider three subsets of the general double trap- 
ping problem: the cases where Da ^ Db ^ Dq, where 
Db — and where Da — Dc — 0- In the first two 
cases it will be necessary, to obtain exact results or exact 
bounds, to limit consideration to the regime pBo ^ Pc- 
For the case Da — Dc = 0, however, an exact result for 
Qa{oo) is possible for all values of the ratio pbo/pc- 
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A. The case Da < Db < Dc 

We first consider a subset of the general A + B ^ B, 
B + C ^ C problem in which the diffusion constants of 
the A, B and C particles are subject to the condition 
Da Db ^ Dc- This allows us to approximate each 
process as an independent "target annihilation problem" 
[31, i.e. in the reaction B+C ^ C we treat the B particles 
as static and the C particles as mobile traps, and in the 
A + B B process we consider the B particles as mobile 
traps for a static A particle. We also subject this problem 
to the condition pso Pc to ensure that the B particles 
remain Poisson distributed even at large times, i.e. that 
no clustering of B particles emerges in regions free of C 
particles. The reason for this condition will become clear 
shortly. 

We begin by considering a single, static B particle at 
the origin. The survival probability, Qi, of the B par- 
ticle with only one C present, starting at some x > 0, 
is the solution of Q subject to the boundary conditions 
Qi(0,t) = 0, Qi(oo,t) = 1 and Qi(x,0) = 1. That is 



Qi{x,t) = erf 



where erf(j/) is the error function. With N traps present, 
where N = pc'L, starting at positions Xi uniformly dis- 
tributed in [0, L], the survival probability of the B parti- 
cle is (following the argument in 0) 



N 



QN{t) = n 



— / da;,, erf 



where we have averaged over the starting positions. 
Rewriting the error function in terms of the complemen- 
tary error function, erf(j/) = 1 — erfc(2/), and using the 
fact that that the C particles are independent, we have 



1 / da; crfc 



PaL 



(3) 



Taking the limits N ^ oo and L — > oo, keeping pc fixed, 
and evaluating the integral in Q, gives 



Q{t) = exp ( — y= pc\jDct 



(4) 



We can perform the same calculation on the opposite 
side of the B particle and, since the results are symmetric 
and independent, get the same result. So the full survival 
probability of the B particle is simply the square of 



For the A + B ^ B process we can use the 'toy model' 
introduced in In this model, the traps B are as- 
sumed to disappear randomly in a manner consistent 
with the required density psit). The model will be an 
exact representation of the double trapping reaction pro- 
vided pbo PC: so that no correlations develop in the 
positions of surviving _B-particles. The time-dependence 
of the A-particle survival probability, QA{t), within this 
model is calculated using a similar argument to that out- 
lined above, but with the traps disappearing from the 
system at each time step with a known probability. In 
our case, we can describe this decay of traps using the sur- 
vival probability of a i? particle, Qsit), as given by jSJ. 
The model gives for the survival probability of the A 
particle 



QAit) = exp -2pBo 



Db 



dr 
7172' 



Substituting for Qb{t) from Eq. © gives the result 



QA{t) = exp 




(6) 



Note that 
as static, 
limit Da 

Dc)- 

For the 
hold, but 
according 
does not 
bound on 



this result was derived treating the A particles 
so it will become asymptotically exact in the 
Db (and we have already assumed Db 

case where pBo ^ Pc and Db Dc still 
Da is arbitrary, the "Pascal principle" [lol |. 
to which the A particle survives longest if it 
move, shows that Eq. © provides an upper 
QA{t) for any value of Da- 



I.e. 



B. The case Db = 

For the B + C ^ C reaction, with the B particles re- 
maining static, we can again use the result in Eq. |SJl for 
the _B-particle survival probability. The A + B ^ B pro- 
cess remains, however, a nontrivial problem, consisting 
of A particles diffusing among static traps each of which 
has a survival probability (t) , as shown in fig. ^ Tak- 
ing once again the limit pbo PCi the survival proba- 
bilities of the different i?-particles can be treated as in- 
dependent, but the problem is still nontrivial. We can, 
however, find a lower bound on QA{t) by considering an 
approach based on the span, i?(i), of a random walk - the 
distance from the point of maximum excursion in one di- 
rection to maximum excursion in the opposite direction, 
up to time t. The average value of this quantity is given 
by (see ^) 



Qsit) = exp 



■PC 



Vd^ 



(5) 



{R{t))=4 



DAt 



(7) 



which is the standard result for the one-dimensional tar- 
get problem. 



We formulate the A+B B process as follows. We write 
the infinitesimal change in the survival probability of a 



FIG. 1: For Db = and pso ^ pc, the A particles diffuse 
among static traps which disappear randomly and indepen- 
dently. 



given A particle, cIQa, in terms of the probability that it 
is trapped by a _B particle in the time interval t t + dt 
in which the span of the random walk increases by dR: 



dQA = -pBodRQ{A,B;t) 

= ~PBQdRQAQ{B\A;t), 



(8) 



where psodR is the probability to find a B particle ini- 
tially in the interval dR (treating, as usual, the initial B- 
particle locations as a Poisson process), Q{A, B; t) is the 
joint probability that both the A particle and the B par- 
ticle have survived up to time t, and Q{B\A; t) is the con- 
ditional probability that the B particle has survived until 
time t given that the A particle has survived. To treat the 
survival probability of the A and the B particles as inde- 
pendent, we must make the assumption that the B parti- 
cles remain Poisson distributed at all times, i.e. that the 
positions of the B particles are not spatially correlated. 
If the B particles survive in clusters and there develop 
regions totally free of B particles, then the probability 
that an A particle has survived will depend on whether 
it is in a i?-free region or in a cluster of B particles. In 
the limit pso ^ PCi the assumption of independence is 
justified, and allows us to write Q{B\A;t) = Qsit), with 
Qsit) given by lO). The solution of ^ is then 

QA{t) = exp (^-pBo J Qb{t) ^ dr 

The span is a stochastic variable so we need to average 
over all realizations of the function i?(r), which analyt- 
ically is not trivial. We can, however, obtain a lower 
bound by using the convexity inequality 

{QA{t)) = ^exp (^-pBQ J Qb{t) ^ dr 

> exp (^-pBQ J Qb{t) ^{R) dr^ ^ 

and by substituting for {R{t)) using {T)) we obtain the 
result 



(QAit)) > exp 




(9) 



L' 







FIG. 2: We consider the subset of the problem where Da = 
Dc = 0. For the A particle to survive, all the B particles 
between the neighbouring C particles must reach a C particle 
first, thus being removed from the system before reaching the 
A particle. 



C. The case Da = Dc = 

We now consider the case where the A and C particles 
remain static and the B particles diffuse among them. 
We can treat this is an extension of the Gambler's Ruin 
problem |^. We need only consider a single A particle 
and the nearest C particle on either side of it, as shown 
in fig. 121 For the probability that the A particle survives, 
we want the probability that all B particles between the 
A particle and the nearest C particles on either side reach 
the C-particles before the A particle. The results (after 
averaging over the distances L and L' in Figure|2I) "^ill be 
the same on each side, and independent so we may solve 
the backward Fokker-Planck equation 101 on one side and 
simply square the result. We solve Eq. Q subject to the 
boundary conditions Q{0,t) = 0, Q{L,t) — 1, where we 
have an A particle at the origin and a C particle at x = L, 
and the B particle starts at x, uniformly distributed in 
[0, L]. The solution is 



Qix,t) 



X 

L 



oo ^ 

E— sin 



n=l 



L 



X exp 



n'^TT'^DBt\ 



i2 



(10) 



Averaging the result over the starting position x gives 



Q{t-L) = Y f dxQ{x,t) 
^ Jo 



1 4 



y2 ~2 



odd n 



n^TT^DBt 
L2 



(11) 



We now consider an arbitrary number of B particles 
in this interval [0,L]. Since the B particles are Poisson 
distributed, the probability to have TV of them initially 
in the interval is 



Pn 



{pbqL) 



N 



-PbqL 
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Then the probabihty that the A survives given N B par- 
ticles initially in [0, L] is, for large t, 



Q{t;L)] 



N=Q 

exp 



— z 1 — e 



where we have kept only the lowest mode since we are 
interested in an asymptotic large-< result. 

Finally, we average over all possible lengths L, 
weighted by the Poisson distribution for the C particle 
positions: 



/>oo 

Ja 



(12) 



We simplify by differentiating with respect to t and eval- 
uate the resulting integral asymptotically for large t using 
the Laplace method [IJ. The result is 



1 



Pbo 
2pc 



valid for p^^DBt 3> 1. To include the contribution from 
the left side we square this result to obtain, asymptoti- 
cally, 



QAit) 



16 



(13) 



Note that this result does not require any condition 
on the ratio pbo/pc- As a check on the result we eval- 
uate ifT^ numerically using Gauss-Legendre two-point 
quadrature. We change variables to u = l/{py + 1), 
where p = pc + Pbo/ "2, to map the infinite range of inte- 
gration onto the finite interval [0,1]. The numerical re- 
sult along with the asymptotic result H13() are compared 
in figs. El and 21 for two values of the ratio pbo/pc, 0-5 
and 2 respectively. For ease of comparison we first write 
the result (|13|) in the form 

g^(t)-g^(oo) ^ exp[-A,(p|oi?,t)V3], (14) 

QA[oo)Xi[p%QDBt)'^''' 

where Ai = 16/%/3^ and A2 = 3iTr/2)'^/^{pc / PBo + 
1/2)^/"^. The result is plotted in log-linear form in figs. El 
andjH 



III. DISCUSSION AND SUMMARY 

The asymptotic result ® for the survival probability 
of an A particle in the case Da Db Dc and the 



Numerical 

- - Asymptotic 



FIG. 3: Convergence of the asymptotic solution 11311 to the 
numerical solution of 11211 with pso/pc = 0.5. The data are 
presented as a log-linear plot in the form suggested by Eq. 
(14). The asymptotic solution has gradient 1. 




FIG. 4: Convergence of the asymptotic solution 11311 to the 
numerical solution of 1121 with pBo/pc = 2. The data are 
presented as a log-linear plot in the form suggested by Eq. 
(14). The asymptotic solution has gradient 1. 



lower bound @ obtained for Db = are both subject 
to the condition pBo ^ pc- Under these conditions both 
results reduce to the same limiting forms as at late 
times: 



. , X , Pbo Da,b 

M 00 = 1 \ -f;—: 

PC V Dc 



(15) 



correct to first order in pbo/pc- Results © and © are, 
however, of a more general form since they indicate the 
nature of the asymptotic time dependence oiQAit)- Un- 
der the condition pBo ^ Pc necessary for these results to 
be valid, the exponentials in © and (O can be expanded 
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to first order in their arguments. 

The result for the condition Dj^ = Dc = is 
still more general since it is valid for any value of the 
ratio pbo/pc- The exact infinite-time result for this 
case, Qa{oo) = (1 + pBof^Pc)"^ , differs from the result 
(1 + PBo/pc)^^ obtained from Eq. |^ under the same 
conditions, although once more the two results reduce to 



the same limiting form, 1 — pbo/pc, to leading order in 
P BO I PC- These results suggest the possibility of a sys- 
tematic expansion in powers of pBo/ pc- 
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